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Abstract 

We consider the generalised Beta function introduced by Chaudhry et al. [J. 
Comp. Appl. Math. 78 (1997) 19-32] defined by 


B(x,y;p ) = 


P -1 (l — t) v ~ x exp 


_4£(1 — t)_ 


dt. 


where > 0 and the parameters x and y are arbitrary complex numbers. The 
asymptotic behaviour of B(x,y;p) is obtained when (i) p large, with x and y fixed, 
(ii) x and p large, (iii) x, y and p large and (iv) either x or y large, with p finite. 
Numerical results are given to illustrate the accuracy of the formulas obtained. 
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1. Introduction 


In [1], Chaudhry et al. introduced a generalised beta function defined by the Euler-type 
integral 1 


B(x,y;p)= /V-gl-tp- 1 

Jo 


exp 


~P 


4f(l - t)_ 


dt, 


( 1 . 1 ) 


where 5f£(p) > 0 and the parameters x and y are arbitrary complex numbers. When 
p = 0, it is clear that when 3fJ(x) > 0 and 9i(y) > 0 the generalised function reduces to 
the well-known beta function B(x,y ) of classical analysis. The justification for defining 
this extension of the beta function is given in [1] and an application of its use in defining 
extensions of the Gauss and confluent hypergeometric functions is discussed in [2]. It is 
evident from the definition in (1.1) that B{x,y,p) satisfies the symmetry property 


B(x,y,p) = B(y,x;p). 


( 1 . 2 ) 


A list of useful properties of B(x, y,p) is detailed by Miller in [4], where it is established 
that B(x,y,p ) may be expanded as an infinite series of Whittaker functions or Laguerre 
polynomials; see (A.l). He also obtained a Mellin-Barnes integral representation for 

x The factor 4 is introduced in the exponential for presentational convenience. 
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B(x,y;p), which we exploit in Section 2, and expressed B(x,x ±rr,p ) and B( 1 ±n, l;p), 
where n is an integer, as finite sums of Whittaker functions. 

Our aim in this note is to derive asymptotic expansions for B(x,y;p ) for large x, y 
and p. We consider (i) \p\ — > oo in | arg p\ < ^tt, with x and y fixed, (ii) x and p large, 
(iii) x , y and p large and (iv) either x or y large, with p finite. The expansion for large 
p is obtained using a Mellin-Barnes integral representation for B(x,y,p), whereas the 
other cases are obtained using the method of steepest descents. 


2. The expansion of B(x,y;p) for large p with x, y finite 


We start with the Mellin-Barnes integral representation given by Miller [4] 




c+ooz 


T(s)T(x + s)T(y + s) 
T{\x+\y+s)T(\x+\y+\+s)^ 


S ds 


( 2 . 1 ) 


valid in ] arg p\ < where c > max{0, — — 3?(y)} so that the integration path 

lies to the right of all the poles of the integrand situated at s = —k, s = —x — k and 
s = — y — k, k = 0,1,2,.... Displacement of the integration path to the left over the 
poles followed by evaluation of the residues (assuming that no two members of the set 
{0, x, y} differ by an integer - thereby avoiding the presence of higher-order poles) yields 
the result that B(x,y,p ) can be expressed as the sum of three 2-^2 (— \p) hyper geometric 
functions; see [4, Eq. (1.6)]. 

Since there are no poles in the half-plane 5?(s) > c it follows that displacement of the 
integration path to the right can produce no algebraic-type asymptotic expansion; see 
[8, §5.4], We can therefore displace the path as far to the right as we please; on such a 
displaced path, which we denote by L, the variable |s| is everywhere large. The ratio of 
gamma functions in the integrand in (2.1) may then be expanded as an inverse factorial 
expansion given by [8, p. 39, Lemma 2.2] 


r(a)T(x + s)T(y + s ) 
T{\x+\y+s)Y{\x+\y+\+s) 


M -1 


X) (“)' ?c f r ( s + Pm(s)T{s - M - i), 

3=0 


where M is a positive integer and pm{s) = 0(1) as |-s| —> 00 in | arg s| < 7 r. The 
coefficients Cj = Cj(x,y) are discussed below where the leading coefficient Co = 1. 

Substitution of the above inverse factorial expansion into the integral (2.1) then pro¬ 
duces 


B(x, y\p) = 2 1 x V 7r 2 


M—l 

E (-L 


27TZ 


T(s - j - ±)p S ds + R m 


where 

Rm = 77- [ pm{s)Y(s- M -\)p~ s ds. 

2m Jl 

The integral may be evaluated by the well-known Cahen-Mellin integral given by (see, 
for example, [8, p. 90]) 


1 

27ri 



(| arg z\ < :j7r, c > — 5R(a)) 



Asymptotics of a generalised Beta function 


3 


to yield 

S M 1 X 

B(x,y-,p) = 2 1 ~ x ~ y irUp~^e~ p ^ {~) 3 Cjp~ 3 + R M >■ 

L 3 =0 3 

A bound for the remainder Rm has been considered in [8, p. 71, Lemma 2.7], from which 
it follows that Rm = 0 (p ^ AI ~2 e~ p ) as |p| —> oo in | arg p\ < ^71. 

Hence we obtain the asymptotic expansion 


B(x,y;p) = 2 1 ~ x ~ y Tr^ p~^e~ p < (-) 

L 3=0 


3 Cjp 3 +0(p M ] 


( 2 . 2 ) 


valid as |p| -A oo in the sector | arg p\ < \tt. The expansion of B(x,y;p) for large p is 
seen to be exponentially small in | arg p\ < ^ 7 r; this is a standard result when there are 
no poles on the right of the path in ( 2 . 1 ) and routine path displacement does not produce 
any useful asymptotic information [ 8 , §5.4]. 

The coefficients Cj for j > 1 can be generated by the algorithm described in [ 8 , §2.2.4], 
It is found that 

ci = \(l + x + y + 2xy -x 2 - y 2 ), 

c 2 = ^(9 + 6(2 + xy)(x + y + xy) - (7 + 4 xy)(x 2 + y 2 ) - 6(x 3 + y 3 ) + x A + y A + 14 xy), 

which are symmetrical in x and y as required by (1.2). A closed-form representation for 
Cj is derived in the appendix, where it is shown that Cj can be expressed in terms of a 
terminating 3 - 72 ( 1 ) hypergeometric function given by 


cj = Cj{x,y) = 


A)j(y + 




■ .3 -74 


-j, \y ~ \x, \y-\x + \ 


\^y + h 


(2.3) 


where (a)j = T(a + j)/r(a) is the Pochhammer symbol. When x 
simpler expression 


Cj(x, x) 


+ \)j 

j! 


y, this reduces to the 


(2.4) 


We remark that the asymptotic expansion of B(x,y,p) for p —> 00 could also have 
been obtained by application of the method of steepest descents, which we shall employ 
in the subsequent sections. See also the appendix for a different approach. 


3. The expansion of B{x,y\p ) for large x and p with y finite 


We consider the expansion of B(x, y. p) for large x and p. with y finite, when it is supposed 
that p = ax, where a > 0 and | arg x\ < \tt. By the symmetry property (1.2), the same 
result will also cover the case of large y and p, with x finite. From (1.1), we have 


where 


P 1 

B(x,y;ax ) = / f dt (|arg x\ < i-7r), 

J 0 

(1 ~t) y 


(3.1) 


iit) = 


- log t, f(t) = 


4t(l-t) ” ' ^ v 7 t 

Saddle points of the exponential factor are given by = 0; that is, at the roots of the 
cubic 


f(l - t) 2 + ±a(l - 2f) = 0. 


(3.2) 
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We label the three saddles to, t\ and 1 2 . All three saddles lie on the real axis with to 
situated in the closed interval [0,1], with t\ > 1 and t -2 < 0. The f-plane is cut along 
(—oo,0]. Paths of steepest descent through the saddles t r (r = 0,1) are given by 

$${e l6 (ip(t) — i/j(t r )} = 0 , 9 = arg x\ 

these paths terminate at t = 0 and t = 1 in the directions \6 — cj)\ < and \n < 6 — <j) < 
§7r, respectively, where 4> = arg t. 

When x > 0, the integration path coincides with the steepest descent path over the 
saddle to', for complex x in the sector | arg x\ < \ix, the steepest descent path through to 
becomes deformed but still terminates at t = 0 and t = 1; see Fig. 1. Application of the 




Figure 1: The steepest descent and ascent paths through the saddles to and t\ (heavy dots) when 
a = 1/3 and (a) 9 = 0 and (b) 6 = 7 t/4 . The arrows indicate the integration path. In (b) the 
steepest ascent paths spiral round t = 0 out to infinity passing onto adjacent Riemann surfaces. 
The saddle t 2 on the branch cut on (— 00 , 0] is not shown. 


saddle-point method then yields the leading behaviour 


B(x, y; ax) ~ 


2 vr 


xi/}"{t 0 ) 


f{t o)e 


-xip(t 0 ) 


2t r 


Xlj}"{to) 


to 1 {^-to) y 1 exp 


—ax 


4fo(l — to) 


(3.3) 


as |x| —> 00 in the sector 
shows that 


arg x\ < t}7t, where some routine algebra combined with (3.2) 


nto) 


1 — 3to + 4tg 

t§(l-to)( 2 t 0 -l)’ 


We remark that the saddle to = to(a) has to be computed for a particular value of the 
parameter a, either directly from (3.2) or as a cubic root. 

The asymptotic expansion of B(x,y,ax ) is given by [7, p. 47] 


D , s 0 —xib(tn) C2nF(?7. + 2 ) 

B(x, y- ax) ~ 2e x ^ to > ^- 2 


(|x| —>• 00 , | arg x\ < bn). 


n=0 


(3.4) 
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The coefficients C n can be obtained by an inversion process and are listed for n < 8 in [3, 
p. 119] and for n < 4 in [9, p. 13]. Alternatively, they can be obtained by an expansion 
process to yield Wojdylo’s formula [10] given by 


C — 1 V/, ( yC2 n +\)j R . 

n ““ 9„(n+l)/2 ^ n ~ k ^ M J ’ 

Za 0 k =0 j =0 j- a 0 


(3.5) 


see also [5, 6]. Here Byj = Bkj(a±,a 2 , ■ ■ ■ ,a,k-j+ i) are the partial ordinary Bell polyno¬ 
mials generated by the recursion 2 


k-j+i 

Bkj — ] O'rBk—rj— i, A/,-0 — 4o, 

r —1 

where 5 mn is the Kronecker symbol, and the coefficients a r and b r appear in the expansions 


~^(t 0 ) = J2 a r (t-to) r+2 , f(t) = J2br(t-t 0 ) r (3.6) 

r=0 r= 0 

valid in a neighbourhood of the saddle t = t 0 . 

In numerical computations we choose a value of the parameter a and compute the 
saddle to from (3.2). With a value of y , Mathematica is used to determine the coefficients 
a r and b r for 0 < r < no- The coefficients C 2 n can then be calculated for 0 < n < no 
from (3.5). We display the computed values of C 2 n for different values of a and y in Table 
1. In Table 2, the values of the absolute relative error in the computation of B(x,y,ax ) 
from (3.4) are presented as a function of the truncation index n when x = 100. 


Table 1: Values of the coefficients GW (to lOdp) for different a and y. 


n 

a = 1, y = 1 

« = 5 , v = 1 

a = §> y = l 

a = 2, y = \ 

0 

+0.2668661228 

+0.1364219142 

+0.2036093538 

+0.3909054941 

1 

+0.0982652355 

+0.2683838462 

+0.0762869817 

-0.0309094064 

2 

-0.0635656655 

-0.1085963949 

-0.0456489054 

-0.0039290992 

3 

+0.0186002666 

+0.0151339630 

+0.0137423943 

+0.0024209801 

4 

-0.0039253710 

-0.0003383888 

-0.0026770977 

-0.0005115807 

5 

+0.0012059654 

+0.0004533741 

+0.0003423270 

+0.0000299402 


4 . The expansion of B{x,y]p) for large x , y and p 


We consider the expansion of B(x,y,p ) for large x , y and p, when it is supposed that 
p = ax and y = bx, where a > 0, b > 0 and | arg x\ < From (1.1), we have 


B(x,yp)= fme-^dt 
Jo 


(| arg x\ < ^7r) 


(4.1) 


2 For example, this generates the values B41 = <24, B42 = a 2 + 2aia3, B43 = 3a 2 a2 and B44 = af. 
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Table 2: Values of the absolute relative error in B(x, y, ax ) when x = 100 for different truncation index. 


n 

a = 1, y = 1 

a =by= 1 

a = b v = t 

a = 2, y = \ 

0 

1.838 x 10 -3 

9.682 x 10 -3 

1.853 x 10 -3 

3.963 x 10 -4 

1 

1.770 x 10" 5 

5.892 x 10" 5 

1.666 x 10" 5 

7.426 x 10" 7 

2 

1.295 x 10~ 7 

2.058 x 10" 7 

1.255 x 10" 7 

1.153 x 10" 8 

3 

9.506 x 10” 10 

1.517 x 10" 10 

8.562 x 10" 10 

8.568 x 10" 11 

4 

1.295 x 10' 11 

9.526 x 10" 12 

5.011 x 10" 12 

2.332 x 10" 13 

5 

3.688 x 10" 12 

1.933 x 10 -13 

5.472 x 10 -14 

6.917 x 10 -15 


where 

m = yTT) - log ‘ - Mog<1 - /(f) = 

Saddle points of the exponential factor are given by the roots of the cubic 


t( 1 - t){ 1 - {b + 1 )t} + \a{ 1 - 2 t) = 0. 


(4.2) 

(4.3) 


Routine examination of this cubic shows that, when a > 0, b > 0, all roots are real, with 
one root greater than 1, one in the interval [0,1] and one negative root. The distribution 
of the saddles is thus similar to that in Section 3, where we continue to label the saddle 
situated in [0,1] by to- The topology of the path of steepest descent through the saddle 
to, given by Q{e ie (i/j(t) — ^(to)} = 0 where 6 = arg x, is also similar to that depicted in 
Fig. 1. 

Accordingly, the expansion of B(x,y,p) when p = ax and y = bx, with a > 0, b > 0, 
is given by 

B(x, bx\ ax) ~ 2 e~ x ^ho) " 2n [ — _t ll (|x| —>■ oo, |argx|<^7r), (4.4) 

71=0 X n+ 2 


where the coefficients C-^n can be determined from (3.5) when the coefficients a r and b r 
in (3.6) are evaluated from the definitions of ip(t) and /(t) in (4.2). 

The leading behaviour is 


B(x, bx\ ax) 


2vr 


xip"(to) 
I 27T 


f{to)e 


-X1p(to) 


as x’ 


V x^"(t 0 ) 0 
oo in the sector | arg x\ < ^vr, where 
1 — 3to + 4 1 


tVil-to) 


bx—1 


r(t 0 ) = 


i - 


bt 0 


tg(l — t 0 )(2t 0 — 1) V 1 to J to(l-to) 2 


exp 


+ 


—ax 


L 4to(1 — to)J 


and to = to (a, b) is the root of (4.3) situated in t € [0,1]. 

We note that when b = 1 we have the result [1, 4] 

B(x,x-,p) = 2 1 - 2x Tr^p^- 1) / 2 e-^W_i x i x (p) 


(4.5) 


in terms of the Whittaker function see (A.l). 
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5. The behaviour of B(x,y,p) for large x and finite y and p 

In this final section, we examine the behaviour of B(x,y 1 p) for large complex x = \x\e ie , 
with 0 < 9 < 7r, when y and p > 0 are finite. The situation when — n < 9 < 0 is 
analogous and, in the case of real y, B(x,y;p ) assumes conjugate values. This case 
has been discussed in [2, Appendix], but is repeated (with minor corrections) here for 
completeness. By the symmetry property (1.2), the same result will also cover the case 
of large y, with x and p finite. 

From (1.1), we have upon interchanging x and y (by virtue of (1-2)) 

B(x,y-p)= f (5.1) 

Jo 

where 

«*>-*(rV e " log(1 - () - /(1) = ru' < 5 - 2 > 

Because p > 0 is a fixed parameter, the integral (5.1) is valid for arbitrary complex values 
of x and y. Saddle points of the exponential factor arise when ijj'it) = 0; that is, when 

t 2 (t — 1) + ae~ ie {1 — 2t) = 0. (5.3) 

We label the three saddles to, ti and ^2 as in Section 3. When 9 = 0 , all three saddles are 
situated on the real axis with to € [0,1] and t\ > 1, t 2 < 1- As 9 increases, the saddles 
to and t 2 rotate about the origin and t\ rotates about the point t = 1. The result of this 
rotation is that, when 9 = 7r, to and t 2 become a complex conjugate pair near the origin 
and t\ is situated in the interval [0,1]; see Fig. 2. 






Figure 2: The steepest descent and ascent paths through the saddles to and t\ (heavy dots) 
when a = 1/3 and (a) d = 0.25t r, (b) 6 = 9 0 = 0.65595tt, (c) d = 0.69tt, (d) 6 = 6 1 = 0.71782tt, 
(e) 6 = 0.807T and (f) 8 = ir. The arrows indicate the integration path. The steepest ascent paths 
spiral round t = 1 out to infinity passing onto adjacent Riemann surfaces. The saddle t 2 is not 
shown. The t-plane is cut along [l,oo). 
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When 0 = 0, the integration path coincides with the steepest descent path passing 
over the saddle to given approximately by 

1 i , \ 

to — a 2 — ^a (x —>■ oo). 

Then, with the estimates 

xifj(to) ^ (px ) 1/2 + | p, ip"(t 0 ) ~ 2 a~^, 
we find by application of the saddle-point method the leading behaviour 


B(x, y;p) 


—y —— 

7T ( p \ 2» 4 


exp 


— (px ) 1 ! 2 — 


;P 


(0 = 0, x -A Too). 


(5.4) 


x \4x / 

When 0 = 7T, we find from (5.3) that the saddle t\ close to the point t = 1 is given by 


ti ~ 1 — a T cr 


( 1*1 


oo 


and 


\x\ ip(t\) ~ |x| T \p — |x| log a, ip"(ti) — a 2 . 

The integration path again coincides with the steepest descent path through t\, and so 
we obtain the behaviour 


(0 = 7T, X —> — OO) 

(5.5) 

The leading terms in (5.4) and (5.5) were given in [2, Appendix], 

A detailed study of the topology of the steepest descent paths 3 through the saddles 
to and t\ when 0 < 0 < it is summarised in Fig. 2 for the particular case a = 4. 
The t-plane is cut along [l,oo) and paths of steepest descent either terminate at t = 0 
(with | arg t\ < \i r), t = 1 (with | arg(l — t)\ < - 771 ) or at infinity. Paths that approach 
infinity spiral round the point t = 1 passing onto adjacent Riemann surfaces. The figures 
reveal that there are two critical values of the phase 0, where the saddles to and t\ 
become connected (via a Stokes phenomenon). We denote these values by 0o = 0o(a) 
and 0i = 0i(a), where a is defined in (5.2). The values of these critical angles are 
tabulated in Table 3 for different a. 

When 0 < 0 < 0o(a), the integration path can be deformed to coincide with the 
steepest descent path passing over to, so that the leading behaviour in (5.4) applies in 
this sector. When 60 (a) < 0 < 0i(a), the integration path is deformed to pass over both 
saddles to and fi, where each steepest descent path spirals out to infinity. Finally, when 
0i (a) < 0 < 7 r, the integration path is deformed to pass over only the saddle t\. 

Based on these considerations and on the approximation of the saddles to — a' 2 — \a!, 
ti ~ 1 + a' - a ,3 , where a! = p/(4x), the leading behaviour of B(x,y;p ) is found to be 


B(x,y\p ) 



exp 


x ~ ~P 


-M 


exp 


-1*1 - 4P 


f J 0 0 < 0 < 0i (a) 

B(x,y,p) ~ I Jo-Ji 61 (a) < 6 < 62 (a) (5.6) 

[ Ji 02(a) < 6 < TT 

3 The saddle T does not enter into our consideration as it plays no role in the asymptotic evaluation 

of B(x,y,p) when 0 < 9 < n. 



Asymptotics of a generalised Beta function 


9 


Table 3: Values of the critical angles 9q, 9\ and 9* as a function of a = p/(4|x|). 


a 

9 0 /n 

9 1 /n 

9*/tt 

0.30 

0.603324 

0.752315 

0.688289 

0.25 

0.536784 

0.798621 

0.681218 

0.20 

0.476795 

0.840611 

0.672858 

0.15 

0.418651 

0.879708 

0.662628 

0.10 

0.358268 

0.916935 

0.649359 

0.05 

0.288029 

0.953688 

0.629820 

0.01 

0.198480 

0.986248 

0.597144 


as \x\ 


and 


oo when 0 < 9 < it (with y and p > 0 finite), where 
Jq : = 


Jl : = 


2vr 


*o 1 1 1 ' exp 


i„-i 
7T ( p \ 2 y 4 


x\V(to) 

P_ 

x \4x 

27r 

2 n 


-P 


4t 0 (l — to)J 


exp 


— (px ) 1 / 2 — 


-,P 


t\ 1 (1 — t\) x 1 exp 


x—1 


-P 


L4ti(l — ti)J 


~ l 


( P Y nix 

1 1 

— e exp 
\4 x ) 

x — —p 


(5.7) 


(5.8) 


with arg ip"(t r ) € [0, 27r], r = 0,1. Inspection of Table 3 shows that as a decreases (that 
is, as |x| increases for fixed p) the angular sector do(a ) < 8 < 6 i(ct), where B(x,y;p ) 
receives a contribution from both saddles, increases. We also show in Table 3 the value 
of 6 = 9*(a) at which ^(^(fo)) = when the saddles are of the same height. We 

have («) < 0*(a) < 0i(ct); then, for 9 < 9* (a) the saddle to is dominant, whereas when 
9 > 9*(a ) the saddle t\ is dominant in the large-|x| limit. 

In Table 4 we present the results of numerical calculations using the asymptotic 
behaviour of B(x, y,p) in (5.6) compared to the values obtained by numerical integration 
of (5.1). The parameter values chosen correspond to a = 0.01 and the saddles to and t\ 
are computed from (5.3), with the leading forms Jo and J\ computed from (5.7) and (5.8). 
It is seen from Table 3 that the exchange of dominance between the two contributory 
saddles arises for 9 ~ 0.607T. 


Appendix: A closed-form expression for the coefficients Cj 


In this appendix we derive a closed-form expression for the coefficients Cj appearing in 
the expansion (2.2). Miller [4, Eq. (2.3a)] has shown that B(x,y,p) can be expressed as 
a convergent series of Whittaker functions in the form 


B(x,y-p) = 2 1 -*-%V y - 1)/2 e-^E ( ^ ^ * g)fc 


W_ k _ i i 


Xp), (A- 1 ) 


k =0 
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Table 4: Values of the asymptotic behaviour of B(x,y,p) in (5.6) with the calculated value when 
|a;| = 50, p = 2 (a = 0.01) and y = ^ for different 9 = arg x. 


0/it 

Asymptotic value 

Calculated value 

0 

0.20 

0.40 

0.50 

0.60 

0.70 

0.80 

1.00 

+5.175 x 10" 06 

-8.210 x 10 -06 + 2.081 x 10 -06 i 
+3.468 x 10“ 05 - 6.934 x 10" 06 z 
+2.647 x 10" 06 - 9.853 x 10" 05 i 
-8.837 x 10" 04 - 3.821 x lO" 03 ? 
-5.944 x 10+ 28 + 1.659 x 10 +28 i 
+2.786 x 10+ 54 + 3.451 x 10 +54 i 
+4.146 x 10+ 77 

+5.187 x 10"° 6 

-8.223 x 10 -06 + 2.096 x lO" 06 ? 
+3.470 x 10"° 5 - 7.020 x 10" 06 i 
+2.402 x 10"° 6 - 9.855 x 10“ 05 i 
-8.781 x 10" 04 - 3.823 x lO" 03 ? 
-5.952 x 10+ 28 + 1.652 x 10 +28 i 
+2.786 x 10+ 54 + 3.459 x 10 +54 i 
+4.154 x 10+ 77 


where W K)/l (x) is the Whittaker function. For p —>• oo with bounded k, we have the 
expansion [7, Eq. (13.19.3)] 


W_ k _i yt i y (p) =p k h e h p 


N 1 n + k) n {y + 5 + k) r 


EH 

n=0 


n\p n 


+ 0 (p 


-AA 


where N is a positive integer. Then we obtain from (A.l) 

B(x,y;p) = 2 1 ~ x - y irh p-h e - p {S(x,y,p) + 0(p~ N )}, 

where 

N ~ 1 ( hy-hx) k (h+hy-hx)k E, 


(A.2) 


S(x, y,p) = 

k =o 


k\p k 


EH 

?r=0 


n\ p n 


= E 

fc=0 


( 2 t/ 2 x)fc( 2 + 2 y 2 ®)fc ^ / M-fc ( 2 4~fc)j-fc(?/+ 


fc! 


EH j - 

j=k 


(j -k)\pi 


+ 0(p 


-JVn 


and we have made the change of summation index n —> j — A;. Use of the fact that 
( —j)k = (—)H/(i — &)!, the above double sum can be written as 

^ 1 (E — E) fc ‘E 1 ( — V r(j+|)r(t/+j+i) 


E 

k =0 
iV-1 


A:! 


E 


j! r(/c+±)r( y +±+*H 




_ (_ Y (2)j(2/ A 2)j y^ (_j)fc( 2 V 2 x ) k ^ 2 ~ >r 2y 2 x ) k 

< —/ A| /> I ml • ^ 


j =0 3’ J ' 1 ' k =0 

= E (H (jMHjk 3F2 


j-p° 


Ho J l 


j\p j 


k\ {\)k{y + H 


j) 2 y 2 ®’ 2~^2 y 2 x 

bv + h 


1 


(A.3) 


upon reversal of the order of summation and identification of the inner sum over k as a 
terminating 3+2 series of unit argument. 

Comparison of (A.2) and (A.3) with the expansion obtained in (2.2) then yields the 
final result 

{\)j{y + \)j ^ I 


-3, h-\x,\ + \y-\x / 


/A A \ 
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